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ABSTRACT 


The  Harpoon  Weapon  System  (HWS)  is  an  all-weather,  over-the- 
horizon,  antiship  missile  deployed  on  U.S.  Navy  vessels  and 
aircraft.  Shipborne  versions  of  the  HWS  are  comprised  of  the 
Harpoon  missile,  missile  launchers,  and  the  Harpoon  Ship  Command- 
Launch  Control  System  (HSCLCS) .  In  order  to  add  the  capability 
to  display  geographic  information  on  the  HSCLCS,  the  Naval 
Research  Laboratory  (NRL)  has  developed  databases  and 
retrieval/display  algorithms  for  both  vector  shoreline  data  and 
gridded  landmask  data.  The  shoreline  database  used  is  a  thinned 
and  compressed  version  of  the  Defense  Mapping  Agency's  (DMA) 

World  Vector  Shoreline  (WVS) .  The  landmask  database  is  a  1- 
minute  grid  of  land/water  values  derived  from  the  WVS.  Extensive 
guality  control  procedures,  including  paper  plots  of  the  entire 
world  and  workstation  display  of  all  shoreline  regions,  has 
ensured  that  the  two  databases  are  in  agreement  and  that  the  data 
are  correct. 
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Interim  Report-Map  Data  Support  for  AM/SWG-1A(V)  Harpoon 
Ship  Command-Launch  Control  System  (HSCLCS) 


BACKGROUND 

The  Harpoon  Weapon  System  (HWS)  is  an  all-weather,  over-the- 
horizon,  antiship  missile  designed  for  use  on  U.S.  Navy  ships, 
submarines,  and  aircraft.  Harpoon  is  now  serving  as  the  Navy's 
basic  antiship  missile  for  Fleet-wide  use,  being  deployed  on 
cruisers,  destroyers,  frigates,  and  patrol  hydrofoils. 

Shipborne  versions  of  the  HWS  are  comprised  of  the  Harpoon 
missile,  missile  launchers,  and  the  Harpoon  Ship  Command-Launch 
Control  System  (HSCLCS) .  The  HSCLCS  configuration  now  in  use, 
AN/SWG-1A(V) ,  has  a  Graphic  Display  System  that  can  display  the 
designated  target  and  other  hostile,  friendly,  or  neutral  ships 
in  the  area  of  interest.  The  HSCLCS  will  automatically  generate 
and  display  an  optimized  flight  path  to  the  target.  The  operator 
can  immediately  launch  the  missile  utilizing  this  flight  plan  or 
he  can  alter  this  missile  flight  plan  in  order  to  conceal  the 
launch  point,  fly  around  other  ships  or  landmasses,  avoid 
defensive  systems,  or  approach  the  target  from  multiple 
directions. 

The  HSCLCS  is  being  upgraded  to  provide  coastline  and 
landmass  display  capability.  The  HSCLCS  currently  performs  the 
functions  of  mission  planning,  launch,  and  control  without 
displaying  landmasses  or  accounting  for  them  during  automatic 
planning.  Effects  of  land  must  be  manually  assessed  and  the 
flight  plan  modified  by  the  operator.  In  the  near  term,  the 
HSCLCS  is  being  upgraded  to  provide  landmass  and  coastline 
display  as  visual  aids  for  engagement  planning.  Future  plans 
call  for  integration  of  the  map  data  into  the  automatic 
engagement  planning  algorithms  to  ensure  land  avoidance. 

The  Graphic  Display  System  provides  a  9.6"  high  x  10.7"  wide 
plasma  display  with  60  pixels/inch  resolution.  The  graphic 
display  is  monochrome  at  range  scales  of  256,  192,  128,  64,  32, 
and  16  kiloyards  radius,  with  a  viewing  area  limited  to  a  four 
inch  radius  circle.  The  32  kiloyard  scale  is  used  for  most  search 
pattern  assessments,  with  the  16  kiloyard  scale  selectable  if 
required. 

The  upgraded  HSCLCS  will  use  an  Intel  80486  microprocessor, 

8  megabytes  of  high-speed  RAM,  and  24  megabytes  of  high-speed, 
nonvolatile  memory.  The  Naval  Research  Laboratory  (NRL)  is 
developing  databases  and  retrieval  algorithms  to  support  display 
of  land  areas  anywhere  in  the  world  under  these  system 
constraints.  Geographic  data  will  be  updated  as  required  by 
vessel  movement.  The  databases  are  intended  to  reside  in  flash 
memory  until  required  for  display  updates.  At  that  time. 


geographic  sections  will  be  moved  from  the  high-speed  memory  into 
RAM  and  then  used  to  update  the  display. 

In  order  to  provide  a  high-resolution  geographic  display 
with  sufficient  update  rates  and  limited  screen  clutter,  the 
system  will  use  two  compressed  databases.  Coastline  vectors  will 
be  provided  by  thinning  and  compressing  the  Defense  Mapping 
Agency's  (DMA's)  World  Vector  Shoreline  (WVS) .  Landmask  data 
will  be  provided  by  a  compressed  gridded  database  of  1 ' 
land/water  values  derived  from  the  WVS  data. 

COASTLINES 

The  WVS  is  a  vector  database  developed  and  supported  by  DMA. 
WVS  data  are  derived  primarily  by  vector  processing  Digital 
Landmass  Blanking  (DLMB)  data,  with  some  data  in  polar  regions 
digitized  from  paper  charts  or  taken  from  DMA's  Digital  Terrain 
Elevation  Data  (DTED) .  DLMB  and  DTED  are  gridded  data  bases  with 
a  post  spacing  of  3  arc-seconds  (approximately  100  meters) .  Data 
density  for  the  WVS  averages  approximately  12  data  points  per 
nautical  mile,  or  the  equivalent  of  a  1:250,000  scale  map.  The 
accuracy  requirement  for  WVS  data  is  that  90%  of  all  identifiable 
shoreline  features  be  located  within  500  m  circular  error  of 
their  true  geographic  positions  with  respect  to  the  World 
Geodetic  System  (WGS)  84  datum. 

COASTLINE  THINNING 

WVS  data  were  found  to  have  a  high  percentage  of  production 
"noise"  and  redundant  data  points.  The  data  were  thinned  by  the 
Naval  Research  Laboratory  using  the  Douglas-Peucker  algorithm 
with  a  thinning  threshold  of  2  arc-seconds  (Landrum,  1989) .  The 
2  arc-second  (approximately  60  m)  thinning  is  essentially  a 
noise-removal  process  and  results  in  a  decrease  in  file  size  of 
about  28%  with  no  significant  decrease  in  data  quality. 

Differences  in  display  of  the  original  DMA  WVS  data  and 
display  of  the  thinned  data  set  on  the  Graphic  Display  System 
will  not  exceed  1  pixel  width  (0.0167  inches)  for  a  16  kiloyard 
display.  Differences  at  the  256  kiloyard  range  scale  are  on  the 
order  of  0.001  inches. 

COASTLINE  COMPRESSION 

After  thinning,  the  data  were  compressed  using  the  Depner- 
Hammack  compression  method  described  in  A  Portable  Method  for 
Compression,  Storage  and  Retrieval  of  High-Resolution  Geographic 
Data  Sets  (Hammack,  1990) .  This  compression  technique  is 
lossless  and  allows  the  storage  of  a  single  geographic  data  point 
in  about  2.5  bytes.  The  resulting  WVS  coastline  database  size  is 
about  18  megabytes. 
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FILE  STRUCTURE 


The  coastline  database  file  contains  two  major  sections. 

The  first  is  a  bitmap  that  indicates  the  1°  cells  that  contain 
coastline.  The  second  section  contains  the  actual  latitude  and 
longitude  coordinates  of  the  coastline  segments.  The  coastline 
file  structure  is  shown  in  Figure  l. 

Algorithm  addressing  is  used  in  this  compression  method  to 
provide  a  unique  record  number  for  every  1°  cell  on  the  surface 
of  the  Earth.  The  1°  cells  are  defined  by  the  latitude  and 
longitude  of  their  lower  left  (southwest)  corner.  The  record 
number  is  computed  by: 

record_number  =  (latitude  degrees  +90)  *  360  +  (longitude  degrees  +  180)  +  1 

where  latitude  degrees  is  the  integer  portion  of  the  latitude 
(south  is  negative)  and  longitude  degrees  is  the  integer  portion 
of  the  longitude  (west  is  negative) .  This  formula  yields  a 
record  number  between  1  and  64,800  for  any  cell  on  the  surface  of 
the  Earth. 

The  majority  of  1°  cells  on  the  surface  of  the  Earth  do  not 
contain  any  coastline  information.  A  64,800-bit  bitmap  (1  bit 
per  1°  cell)  is  used  to  indicate  which  cells  contain  coastline. 
This  bitmap  is  stored  in  the  records  2  through  23.  The  records 
corresponding  to  cells  that  do  not  contain  coastline  are  used  to 
store  overflow  records  from  other  cells.  In  this  way  storage  is 
not  wasted  by  empty  records. 

The  use  of  a  1°  bitmap  provides  a  method  for  quickly 
determining  whether  coastline  data  are  available  in  a  cell. 

After  the  bitmap  is  checked  to  see  if  any  data  are  available  in  a 

1°  cell,  the  record  number  is  used  to  directly  access  the  data 
stored  in  the  second  section  of  the  file. 
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Logical  Record 
L*  Bytes: 


Physical  Record 


P*  Bytes 


WWWWWWWW 


First  Logical  Record  (next  empty 
record  pointer,  logical  record  size, 
version  number) 


Cell  Bit  Map  (length  is 
648(X)/(L*x)  + 1  logical  records 

Map  Data  (initial  and  overflow 
records) 

Overflow  Records  Past  Cell  64800 


Figure  1.  General  Direct-Access  File  Structure. 
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In  the  second  section  of  the  file,  coastline  segments  are 
stored  as  lines  that  begin  and  end  within  a  single  cell. 

Coastline  segments  which  would  traverse  more  than  one  cell  are 
broken  at  the  cell  boundaries  during  the  compression  process  and 
treated  as  a  series  of  smaller  segments.  Each  of  these  segments 
is  contained  completely  within  one  cell,  with  the  exception  of 
one  overlap  point  at  the  cell  boundary.  This  overlap  point  is 
required  to  provide  a  smooth  connection  to  the  bordering  cell  and 
eliminate  gaps  at  the  cell  boundaries. 

The  computed  record  number  is  used  as  a  pointer  to  the  3072- 
byte  record  containing  the  first  of  the  coastline  coordinates  for 

the  given  1°  cell.  Multiple  records  are  chained  together  to 
provide  sufficient  storage  space  when  required.  If  there  are 
more  data  than  can  be  stored  in  a  single  3 07 2 -byte  record,  an 
"overflow  pointer"  gives  the  address  of  the  next  record  in  the 
chain. 

Within  an  individual  cell,  coastline  segments  are  stored 
using  a  segment  header  and  delta  records.  Each  coastline  segment 
has  a  single  segment  header  and  multiple  delta  records.  The 
segment  header  defines  three  attributes  of  the  line  segment:  (1) 
the  beginning  latitude  and  longitude  of  the  segment,  (2)  the 
number  of  points  in  the  segment,  and  (3)  a  "rank"  or  identifier 
for  the  type  of  segment.  For  this  application  all  coastline 
segments  are  of  the  same  type,  so  the  "rank"  is  not  used.  The 
structure  of  the  segment  header  is  shown  in  Figure  2. 

The  segment  header  is  6  bytes  in  length.  The  first  bit  of 
byte  1  contains  the  continuation  flag.  If  set,  this  indicates 
that  this  record  is  a  continuation  of  a  previous  record  in  the 
same  cell.  The  remaining  seven  bits  of  BYTE  1  contain  the  most 
significant  seven  bits  of  the  segment  length  (number  of  points) 
and  the  first  four  bits  of  BYTE  2  contain  the  least  significant 
four  bits  of  the  segment  length. 

The  degrees  portion  of  the  starting  latitude  and  longitude 
is  implicit  in  the  record  number  and  need  not  be  stored.  The 
minutes  portion  of  the  starting  latitude  and  longitude  is  biased 
by  128,  so  that  negative  values  may  be  stored,  and  stored  in 
BYTES  5  and  6.  It  is  necessary  to  store  negative  values  in  order 
to  accommodate  the  overlap  points  needed  when  a  coastline  segment 
crosses  a  cell  boundary.  (Negative  values  would  be  required  if 
the  line  segment  exits  the  west  or  south  side  of  the  cell.)  The 
seconds  portion  of  the  latitude  is  split  between  BYTES  2  and  3  as 
shown  in  Figure  1,  and  the  seconds  portion  of  the  longitude  is 
stored  in  a  portion  of  BYTE  3. 
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SEGMENT  HEADER: 


1 

2 

3 

4 

5 

6 

1 

1 

1 

RANK 

LATITUDE 

LONGITUDE 

1 

1 

1 

MINUTES  1 

MINUTES 

BYTE  1  of  the  SEGMENT  HEADER; 
Bit: 


1  2 

3 

4 

5 

6 

7 

00 

Cont. 

Flag 

1 

^  ! 

MOST  SIGNinCANT  7 

BITS  OF  TH 

i  SnGMFNT 

LnNGTlI 

BYTE  2  of  the  SEGMENT  HEADER: 


Bit: 

1  2  3  4  5  6  7  8 


MOyr  SIGNIFICANT  4  BITS  OF  TOE 

I  lATITUDE  I  SECONDS  | 

BYTE  3  of  the  SEGMENT  HEADER; 


Bit: 

1  2 

3 

4 

5  6 

7 

8 

1 

~i  1  r 

2  BITS  OF 

LONGITUDE  SECONDS 

J _ 

_J _ ^ _ L 

1 _ 

Figure  2.  Segment  Header  Data  Structure. 
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DELTA  RECORDS 


Following  the  segment  header  is  a  series  of  delta  records. 
Each  delta  record  consists  of  two  bytes;  a  latitude  delta  and  a 
longitude  delta.  Each  delta  is  the  offset  in  seconds  from  the 
previous  point.  These  delta  records  are  also  biased  by  128  to 
accommodate  negative  values. 

The  last  delta  value  for  most  line  segments  will  not 

correspond  to  the  last  point  in  the  1°  cell.  These  last  values 
are  normally  overlap  points,  and  will  also  be  the  first  point  of 
a  new  segment  in  an  adjacent  1°  cell. 

Details  of  the  compression  method  and  data  formats  are 
included  in  A  Portable  Method  for  Compression,  Storage,  and 
Retrieval  of  High-Resolution  Geographic  Data  Sets  (Depner  & 
Hammack,  1988) . 

Retrieval  and  display  software  has  been  written  in  FORTRAN 
77  and  ANSI  C  and  is  included  as  Appendices  A  and  B, 
respectively. 

LAMDMASK 

The  landmask  database  consists  of  a  simple  1-minute  grid  of 
land  and  water  values  derived  from  the  WVS  database.  The  Naval 
Ocean  Systems  Center  (NOSC)  used  the  WVS  database  to  develop  a  3 
arc-second  grid  of  country  codes  (Houghton,  1989) .  This  database 
was  thinned  to  provide  a  1-minute  grid  of  land/water  values,  then 
compressed  to  about  720  kilobytes  from  the  original  49  megabyte 
NOSC  file. 

LANDMASK  THINNING 

The  first  step  in  developing  the  1'  landmask  file  was  to 
thin  the  original  3  arc-second  data  to  1  minute  and  remove 
country  code  information.  The  NOSC  data  were  run-length  coded  by 
country.  The  NRL  thinning  algorithm  sampled  the  NOSC  database  at 
1-minute  increments  and  translated  country  codes  into  a  single 
••land”  code.  This  thinning  resulted  in  a  30  megabyte  file  of 
land  and  water  values,  with  each  bit  corresponding  to  a  1-minute 
cell.  These  data  were  overlayed  with  WVS  data  on  paper  using  a 
CalComp  electrostatic  color  plotter  in  order  to  detect  any  bad 
data  values  in  the  file.  The  file  was  then  manually  edited  on  a 
Sun  SPARCstation  2  in  an  X-window  environment  to  correct  bad  data 
points.  The  final  data  file  was  then  again  plotted  on  the 
CalComp  plotter  as  a  quality  control  check. 
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LANONASK  COMPRESSION 


After  final  edit  and  quality  control  of  the  uncompressed 
data,  the  database  was  compressed  using  a  technique  that  combines 

bitmaps  and  run-length  coding.  Nearly  90%  of  the  1°  cells 
covering  the  Earth's  surface  do  not  contain  shoreline.  The  3600 

1-minute  cells  within  each  of  these  1°  cells  may  be  represented 
by  a  single  land/water  value,  resulting  in  a  significant  decrease 
in  storage  requirements.  One  technique  for  representing  these 
grids  is  to  create  bitmaps  for  larger  geographical  areas  so  that 
many  1-minute  cells  may  be  represented  by  a  single  value. 

The  file  structure  used  here  utilizes  both  10°  bitmaps  and 

1°  bitmaps.  If  a  given  10°  area  is  entirely  land  or  water,  then 
this  information  is  obtainable  from  the  bitmap  and  no  further 

access  of  the  database  is  necessary.  If  a  10°  square  contains 
both  land  and  water,  it  is  then  necessary  to  retrieve  the 

appropriate  1°  bitmap.  However,  if  the  1°  square  in  question 
is  entirely  land  or  water,  no  further  access  is  required.  It  is 

only  when  a  1°  square  contains  both  land  and  water  that  access 
to  the  run-length  coded  information  is  required.  Using  this 
method,  the  majority  of  database  inquiries  require  two  or  fewer 
reads,  even  if  geographic  points  are  accessed  randomly. 

The  logical  record  size  used  for  the  landmask  database  is  32 
bytes.  The  first  four  bytes  of  the  first  record  contain  the 
number  of  1'  land  cells  in  the  database.  The  fifth  byte  contains 
the  database  revision  level.  The  Harpoon  database  is  identified 
by  revision  level  2.  Level  1  data  are  based  on  the  Central 
Intelligence  Agency  (CIA)  World  Data  Bank  II  coastline,  and  Level 
0  is  not  used.  The  following  bytes  provide  a  textual  description 
of  the  database;  in  this  case,  "NRL  1'  WVS  Landmask." 

The  following  81  logical  records  (2592  bytes)  in  the  file 

contain  the  10°  bitmap.  Each  four-byte  integer  contains  either 
an  address,  if  both  land  and  water  are  present  in  the  square,  or 
a  flag  value  indicating  that  the  square  is  either  all  land  or  all 
water.  In  keeping  with  long-standing  military  tradition,  the 
flag  is  set  to  "1  if  by  land,  2  if  by  sea."  The  addresses 

contained  in  this  section  point  to  a  1°  square  map  record. 

Each  1°  bitmap  contains  a  two-bit  flag  for  each  1°  square 
within  its  parent  10°  square.  This  flag  indicates  whether  the  1° 
square  is  all  water  (2),  all  land  (1),  or  both  (0).  Immediately 

following  this  1°  bitmap  is  the  "initial"  run-length  coded  data 
record  for  the  first  1°  square  that  contains  both  land  and 

water.  The  address  of  any  initial  1°  square  record  can  be 
computed  from  the  number  of  "0"  records  preceding  it  (address  = 
number  of  preceding  "0"  records  *  32) . 
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If  the  run-length  coded  data  for  an  entire  1°  square  do  not 
fit  in  one  logical  record  (32  bytes),  then  the  last  two  bytes  of 
the  "initial"  record  contain  a  pointer  to  an  "overflow"  record. 
Similarly,  overflow  records  may  also  contain  pointers  to  other 
overflow  records.  All  overflow  records  are  stored  after  the 
block  of  "initial"  records  in  the  file  (see  Figure  3) . 

The  1°  run-length  coded  records  were  built  by  dividing  the 
square  into  60,  1-minute  latitude  bands,  then  placing  these 
latitude  bands  end-to-end  (starting  from  the  southern  edge)  to 
produce  a  3600  point  array.  This  array  is  then  compressed  using 
run-length  coding.  Each  logical  piece  (sub-record)  of  the  run- 
length  record  may  be  either  one  or  two  bytes  in  length,  depending 
on  the  size  of  the  run  stored.  This  sub-record  consists  of  a 
flag  indicating  the  number  of  bits  used  to  store  the  run-length 
(either  6  or  14  bits) ,  a  flag  that  indicates  whether  the  run  is 
land  or  water,  and  the  run  length.  The  first  flag  is  1  bit  in 
length  and  has  a  value  of  0  if  the  run  length  is  stored  in  6 
bits,  or  a  value  of  1  if  the  run  length  is  stored  in  14  bits. 

The  second  flag  is  1  bit  in  length,  and  is  set  to  0  if  the  run  is 
water  or  1  if  the  run  is  land.  The  following  6  or  14  bits 
contain  a  count  of  the  number  of  points  in  the  run.  If  there  are 
too  many  data  to  fit  into  one  logical  record  (32  bytes) ,  then  the 
last  2  bytes  of  the  record  will  contain  a  pointer  to  the  next 
overflow  record.  This  pointer  is  stored  as  an  offset  from  the 
current  position.  Overflow  records  are  placed  immediately  after 

the  initial  1°  run-length  coded  records  in  the  file,  but  prior  to 
the  next  1°  bitmap  record  (see  Figure  3). 

RETRIEVAL  AND  DISPLAY  SOFTWARE 

Retrieval  and  display  software  was  developed  for  both  the 
WVS  and  landmask  data  files.  Retrieval  code  was  written  both  in 
ANSI  C  and  in  FORTRAN  77,  and  sample  display  code  was  provided 
for  the  IBM  PC/MS-DOS  environment.  This  software  allows  rapid 
access  to  any  1°  cell  on  the  surface  of  the  Earth. 

SUMMARY 

In  order  to  add  the  capability  to  display  geographic 
information  on  the  HSCLCS,  NRL  has  developed  databases  and 
retrieval/display  algorithms  for  both  vector  shoreline  data  and 
gridded  landmask  data.  The  shoreline  data  used  is  a  thinned  and 
compressed  version  of  DMA's  WVS.  The  landmask  data  is  a  1-minute 
grid  of  land/watei  values  derived  from  the  WVS.  Extensive 
quality  control  procedures,  including  paper  plots  of  the  entire 
world  and  workstation  display  of  all  shoreline  regions,  have 
ensured  that  the  two  databases  are  in  agreement  and  that  the  data 
is  correct. 
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Figure  3.  Landmask  File  Structure. 
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FUNCTION  MASKOl  :  RETURNS  A  LOGICAL  FLAG  INDICATING  THAT  THE 
ONE  MINUTE  GRID  NEAREST  TO  THE  LATITUDE  AND  LONGITUDE  POSITION 
PASSED  TO  THE  SUBROUTINE  IS  EITHER  ON  LAND  OR  WATER.  THE  POINT 
MAY  BE  DEFINED  AS  DEGREES  DECIMAL  OR  DEGREES  AND  MINUTES  DECIMAL. 
THE  FLAG  RETURNED  IS  SET  TO  .TRUE.  FOR  LAND  OR  .FALSE.  FOR  WATER. 
THIS  ROUTINE  CALLS  MASK60  TO  GET  THE  DEGREE  SQUARE  MASK. 

THIS  PROGRAM  USES  A  DATA  BASE  DESIGNED  BY: 

JAN  C.  DEPNER  JAMES  A.  HAMMACK 

NAVAL  OCEANOGRAPHIC  OFFICE  NAVAL  OCEAN  R&D  ACTIVITY 

SYSTEMS  TECHNOLOGY  BRANCH  CODE  117 

STENNIS  SPACE  CENTER,  MS  STENNIS  SPACE  CENTER,  MS 

39522-5001  39529- 

AUTHOR  :  JAN  C.  DEPNER,  04/03/89 

ARGUMENTS  : 

LATD  -  (R)  LATITUDE  DEGREES  OF  THE  POINT  (SOUTH  NEGATIVE) 

LATM  -  (R)  LATITUDE  MINUTES  OF  THE  POINT  (UNSIGNED  UNLESS 

LATD  IS  0.0) 

LOND  -  (R)  LONGITUDE  DEGREES  OF  THE  POINT  (WEST  NEGATIVE  OR 

0  TO  360) 

LONM  -  (R)  LONGITUDE  MINUTES  OF  THE  POINT  (UNSIGNED  UNLESS 

LOND  IS  0.0) 

MASKOl  -  (L)  POINT  FLAG  RETURNED 

VARIABLES  : 

MASK  -  (C)  ARRAY  CONTAINING  THE  ONE  MINUTE  GRID  OP  LAND 

MASK  POINTS  FOR  THE  ONE  DEGREE  CELL  THAT  CONTAINS 
THE  LATITUDE  AND  LONGITUDE  POINT  REQUESTED 
LATP  -  (I)  LATITUDE  INDEX  INTO  MASK  ARRAY 

LONP  -  (I)  LONGITUDE  INDEX  INTO  MASK  ARRAY 

LAT  -  (R)  LATITUDE  IN  DECIMAL  DEGREES 

LON  -  (R)  LONGITUDE  IN  DECIMAL  DEGREES 

LOGICAL  FUNCTION  MASKOl (LATD, LATM, LOND, LONM) 

IMPLICIT  INTEGER  (A-Z) 

REAL  LATD, LATM, LOND, LONM, LAT, LON 
CHARACTER*!  MASK(0; 59, 0: 59 ) , ALFLAG, FLAG60 
DATA  PLAT  /-!/  ,  PLON  /-!/  ,  FLAG60  /’M’/ 

CONVERT  LAT  AND  LON  TO  DEGREES  DECIMAL,  ADJUST  TO  NEAREST  MINUTE 

LAT  =  LATD+(LATM/60.0)*SIGN( 1.0, LATD) +90.0 
LON  =  LOND+(LONM/60.0)*SIGN(1.0,LOND) 

IF  (LON. LT. 0.0)  LON  =  LON  +  360.0 
LAT  =  LAT+. 008333333333333 
LON  =  LON+. 008333333333333 

IF  THIS  IS  IN  THE  SAME  ONE  DEGREE  SQUARE  AS  THE  LAST  ACCESS, 

DO  NOT  READ  A  NEW  MASK 

DLAT  =  INT(LAT) 

DLON  =  INT(LON) 

GET  THE  ONE  DEGREE  MASK 

IF  ( DLON. NE. PLON. OR. DLAT. NE. PLAT)  CALL  MASK60( MASK, LAT-90.0, LON, 
&FLAG60 ) 


A-2 


ooo  o o o o o o o o o o o o o o o o o o o o o o o o o o  o o o  ooo  oooo 


IF  THE  ALL  LAND  OR  WATER  FLAG  IS  NOT  SET,  GET  THE  SINGLE  POINT 
FROM  THE  MASK 

IF  (FLAG60.EQ.’  ’)  THEN 
LATP  =  AMOD(LAT, 1.0) *60.0 
LONP  =  AMOD(LON,1.0)*60.0 
ALFLAG  =  MASK (LONP, LATP) 

ELSE 

ALFLAG  =  FLAG60 
END  IF 

SET  THE  LOGICAL  VALUE  BASED  ON  THE  SETTING  OF  ALFLAG 

IF  (ALFLAG. EQ. 'L* )  THEN 
MASKOl  =  .TRUE. 

ELSE 

MASKOl  =  .FALSE. 

END  IF 

SAVE  THE  PREVIOUS  LAT  AND  LON  DEGREES 

PLAT  =  DLAT 
PLON  =  DLON 
RETURN 
END 

SUBROUTINE  MASK60  :  GIVEN  A  REAL  LATITUDE  AND  LONGITUDE  IN  DEGREES 
(WEST  LONGITUDE  NEGATIVE  OR  0  TO  360  LONGITUDE  AND  SOUTH  LATITUDE 
NEGATIVE)  THIS  SUBROUTINE  RETURNS  A  CHARACTER  ARRAY  CONTAINING  A 
ONE  MINUTE  LAND  MASK  FOR  THE  ONE  DEGREE  CELL  THAT  THE  LATITUDE 
AND  LONGITUDE  POINT  IS  IN.  THE  CHARACTER  ARRAY  CONTAINS  3600 
CHARACTERS  AND  IS  SET  UP  AS  'CHARACTER*!  MASK(0: 59 , 0: 59 ) ’ .  EACH 
CHARACTER  IS  SET  TO  EITHER  'L'  OR  'W* .  IF  ALL  THE  POINTS  ARE 
LAND  OR  ALL  ARE  WATER,  THE  1  CHARACTER  FLAG  'ALFLAG'  WILL  BE  SET 
TO  THE  CORRESPONDING  LETTER.  IF  BOTH  LAND  AND  WATER  POINTS  ARE 
PRESENT,  'ALFLAG'  WILL  BE  SET  TO  '  ' .  THE  MASK  IS  SET  UP  WITH 
POINT  'MASK{0,0)'  BEING  THE  SOUTHWEST  CORNER  OF  THE  CELL  AND  POINT 
'MASK(59,59) '  THE  NORTHEAST  CORNER.  THE  FIRST  SUBSCRIPT  REPRESENTS 
THE  LONGITUDE  MINUTES,  THE  SECOND  IS  LATITUDE  MINUTES  (WITH  THE 
COORDINATE  SYSTEM  SET  UP  AS  0  -  360  LONGITUDE  AND  0  -  180  LATITUDE). 

THIS  PROGRAM  USES  A  DATA  BASE  DESIGNED  BY: 

JAN  C.  DEPNER 
NAVAL  OCEANOGRAPHIC  OFFICE 
SYSTEMS  TECHNOLOGY  BRANCH 
STENNIS  SPACE  CENTER,  MS 
39522-5001 

AUTHOR  :  JAN  C.  DEPNER,  04/03/89 

SUBROUTINE  MASK60 ( MASK , LAT , LON , ALFLAG ) 

IMPLICIT  INTEGER  (A-Z) 

REAL  LAT, LON 
PARAMETER  ( PHYREC=1024 ) 

CHARACTER*!  MASK( 0 ; 59, 0: 59 ), ALFLAG 
COMMON  /ADDRS/  LPERP,REC00,REC01 
DATA  FIRST  /!/  ,  PLAT  /-!/  ,  PLON  /-!/ 

IF  FIRST  TIME  THROUGH,  OPEN  FILE  AND  SET  LPERP 

IF  (FIRST.EQ.l)  THEN 


JAMES  A.  HAMMACK 

NAVAL  OCEAN  R&D  ACTIVITY 

CODE  117 

STENNIS  SPACE  CENTER,  MS 
39529- 
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RECORD  LENGTH  FOR  FILE  IS  PHYREC  BYTES  (CHECK  YOUR  MACHINE  FOR 
RECORD  LENGTH  IN  BYTES  OR  WORDS) 

OPEN ( 10, FILE=' data/ landmask. wvB • ,FORM=' UNFORMATTED' , 

&  ACCESS= ■ DIRECT ' ,STATUS=' OLD' ,RECL=PHYREC) 

LPERP  =  PHYREC/32 
FIRST  =  0 
END  IF 

IF  LAT  >=  90.0  SET  TO  ALL  WATER  AND  RETURN 

IF  (LAT. GE. 90.0)  THEN 
ALFLAG  =  'W' 

CALL  ALLMSK( ALFLAG, MASK) 

RETURN 
END  IF 

CALCULATE  SOUTHWEST  CORNER  OF  CELL 

CLAT  =  INT(LAT+90.0) 

IF  (LON.LT.O)  THEN 

CLON  =  1NT(LON+360.0) 

ELSE 

CLON  =  INT(LON) 

END  IF 

IF  THE  LAT  AND  LON  ARE  NOT  IN  THE  SAME  ONE  DEGREE  SQUARE  AS  THE 
PREVIOUS  LAT  AND  LON,  OR  ALFLAG  =  'M'  (FIRST  CALL  FROM  MASKOl), 
PROCESS  THE  CELL 

IF  (CLON.NE. PLON. OR. CLAT. NE. PLAT. OR. ALFLAG. EQ. 'M' )  THEN 
ALFLAG  =  '  ' 

READ  THE  TEN  DEGREE  SQUARE  MAP  AND  CHECK  FOR  ALL  LAND  OR  WATER 

CALL  RDTEN ( CLAT , CLON , ALFLAG ) 

IF  (ALFLAG. NE.'  •)THEN 

CALL  ALLMSK (ALFLAG, MASK) 

ELSE 

READ  THE  ONE  DEGREE  CELL  MAP  AND  CHECK  FOR  ALL  LAND  OR  WATER 

CALL  RDONE( CLAT, CLON, ALFLAG) 

IF  (ALFLAG. NE.'  ')THEN 
CALL  ALLMSK (ALFLAG, MASK) 

ELSE 

READ  THE  RUN  LENGTH  DATA  FOR  THE  ONE  DEGREE  CELL  AND  BUILD 
THE  MASK 

CALL  BLDMSK(MASK) 

END  IF 
END  IF 

SAVE  THE  LAT  AND  LON  DEGREES 

PLAT  =  CLAT 
PLON  =  CLON 
END  IF 
RETURN 
END 

SUBROUTINE  RDTEN  :  READS  THE  TEN  DEGREE  SQUARE  MAP  AT  THE 
BEGINNING  OF  THE  DATA  BASE.  RETURNS  A  FLAG  INDICATING  IF  THE 
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C  ENTIRE  SQUARE  WAS  EITHER  LAND  OR  WATER.  IF  NOT,  THE  VARIABLE 

C  'RECOl'  IS  SET  TO  THE  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  CELL  MAP. 

C  THE  FLAG  RETURNED  IS  SET  TO  'L'  FOR  LAND  OR  'W'  FOR  WATER. 

C 

C  AUTHOR  :  JAN  C.  DEPNER,  04/03/89 

C 

C  ARGUMENTS  : 

C 

C  LAT 

C  LON 

C  TENFLG 

C 

C  VARIABLES  : 

C 

C  PHYREC 

C 

C  IRECIO 

C 

C  CTYPE 

C  LPERP 

C  RECOO 

C 

C  RECOl 

C  PADDR 

C 

SUBROUTINE  RDTEN( LAT, LON, TENFLG) 

IMPLICIT  INTEGER  (A-Z) 

PARAMETER  ( PHYREC=1024 ) 

CHARACTER*!  IRECIO ( 0: PHYREC-1 ) ,TENFLG, CTYPE ( 2 ) 

COMMON  /ADDRS/  LPERP, RECOO, RECOl 
DATA  PADDR  /-!/  ,  CTYPE  /'L^’W'/  ,  PWORD  /-!/ 

COMPUTE  THE  WORD  (4  BYTE)  POSITION  WITHIN  THE  648  WORD  TEN 
DEGREE  SQUARE  MAP,  MAP  STARTS  IN  THE  2ND  LOGICAL  RECORD 

WORD  *  (LAT/10)*36+LON/10 

IF  THE  WORD  POSITION  HAS  CHANGED  GET  THE  NEW  ONE  DEGREE  CELL 
LOGICAL  ADDRESS 

IF  (WORD. NE. PWORD)  THEN 

COMPUTE  THE  LOGICAL  AND  PHYSICAL  ADDRESS  OF  THE  ONE  DEGREE 
CELL  POINTER  WITHIN  THE  TEN  DEGREE  MAP 

LREC  =  (WORD*4)/32+l 
ADDR  =  INT(LREC/LPERP+1) 

IF  THE  PHYSICAL  ADDRESS  HAS  CHANGED,  READ  A  NEW  RECORD 

IF  ( ADDR. NE. PADDR)  THEN 
READ ( 10 , REC=AD0R ) IRECIO 
PADDR  =  ADDR 
END  IF 

COMPUTE  THE  BYTE  POSITION  OF  THE  ONE  DEGREE  CELL  LOGICAL 
ADDRESS  WITHIN  THE  PHYSICAL  RECORD 

PNT  =  MOD( (WORD+8)* 4, PHYREC) 

BUILD  THE  LOGICAL  ADDRESS  FROM  THE  4  BYTES  AT  'PNT' 

RECOl  =  ICHAR(IREC10(PNT) )*16777216+ICHAR( IRECIO (PNT+1) ) *65536+ 
&  ICHAR( IRECIO (PNT+2) ) *256+ICHAR ( IRECIO ( PNT+3 ) ) 


(R)  LATITUDE  DEGREES  OF  THE  POINT  (0  -  180) 
(R)  LONGITUDE  DEGREES  OF  THE  POINT  (0  -  360) 

( C )  ALL  LAND /WATER  FLAG  RETURNED  ( ' L ’  OR  ' W ' ) 


(I)  CONSTANT  SET  TO  NUMBER  OF  BYTES  IN  A  PHYSICAL 
RECORD  IN  THE  DATA  BASE. 

(C)  ARRAY  CONTAINING  AN  INPUT  RECORD  FROM  THE 
DATA  BASE. 

(C)  ARRAY  CONTAINING  THE  LAND  AND  WATER  FLAGS. 

(I)  LOGICAL  RECORDS  PER  PHYSICAL  RECORD. 

(I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  DATA  RECORD 
(NOT  USED  IN  THIS  ROUTINE). 

(I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  CELL  MAP. 
(I)  PREVIOUS  PHYSICAL  RECORD  ADDRESS. 
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PWORD  =  WORD 
END  IF 

IF  THE  TEN  DEGREE  SQUARE  IS  ALL  LAND  OR  WATER,  SET  THE  FLAG 

IF  (REC01.lt. 3)  THEN 

TENFLG  =  CTYPE{REC01) 

ELSE 

TENFLG  =  ’  ' 

END  IF 
RETURN 
END 

SUBROUTINE  RDONE  :  READS  THE  ONE  DEGREE  CELL  MAP.  RETURNS  A 
FLAG  THAT  INDICATES  IF  THE  ENTIRE  SQUARE  IS  EITHER  LAND  OR 
WATER.  IF  NOT,  THE  VARIABLE  'RECOO'  IS  SET  TO  THE  LOGICAL 
ADDRESS  OF  THE  ONE  DEGREE  DATA  RECORD.  THE  FLAG  RETURNED  IS 
SET  TO  'L'  FOR  LAND  OR  'W'  FOR  WATER. 

AUTHOR  ;  JAN  C.  DEPNER,  04/03/89 

ARGUMENTS  : 

LAT  -  (R)  LATITUDE  DEGREES  OF  THE  POINT  (0  -  180) 

LON  -  (R)  LONGITUDE  DEGREES  OF  THE  POINT  (0  -  360) 

ONEFLG  -  (C)  ALL  LAND /WATER  FLAG  RETURNED  ('L'  OR  'W') 

VARIABLES  : 

PHYREC  -  (I)  CONSTANT  SET  TO  NUMBER  OF  BYTES  IN  A  PHYSICAL 

RECORD  IN  THE  DATA  BASE. 

CELMSK  -  (I)  ARRAY  CONTAINING  THE  LOGICAL  ADDRESSES  OF 

THE  ONE  DEGREE  CELL  DATA  RECORDS. 

IRECOl  -  (C)  ARRAY  CONTAINING  AN  INPUT  RECORD  FROM  THE 

DATA  BASE. 

CTYPE  -  (C)  ARRAY  CONTAINING  THE  LAND  AND  WATER  FLAGS. 

LPERP  -  (I)  LOGICAL  RECORDS  PER  PHYSICAL  RECORD. 

RECOO  -  (I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  DATA  RECORD 

(NOT  USED  IN  THIS  ROUTINE). 

RECOl  -  (I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  CELL  MAP. 

PADDR  -  (I)  PREVIOUS  PHYSICAL  RECORD  ADDRESS. 

SUBROUTINE  RDONE (LAT, LON, ONEFLG) 

IMPLICIT  INTEGER  (A-Z) 

PARAMETER  (PHYREC=1024) 

DIMENSION  CELMSK(0:9,0:9) 

CHARACTER*!  IRECOl (0: PHYREC-1 ) ,ONEFLG,CTYPE (2 ) 

COMMON  /ADDRS/  LPERP, RECOO, RECOl 

DATA  PADDR  /-!/  ,  PRECOl  /-!/  ,  CTYPE  /’L’,'W'/ 

COMPUTE  THE  CELL  POSITION  WITHIN  THE  PHYSICAL  RECORD  AND  THE 
PHYSICAL  RECORD  ADDRESS 

CELL  =  MOD (RECOl, LPERP) *32 
ADDR  =  INT(REC01/LPERP+1) 

IF  LOGICAL  RECORD  HAS  CHANGED,  BUILD  A  NEW  CELL  MASK 

IF  (RECOl.NE. PRECOl)  THEN 

IF  THE  ADDRESS  HAS  CHANGED,  READ  A  NEW  RECORD 

IF  ( ADDR. NE. PADDR)  THEN 
READ ( 1 0 , REC= ADDR ) I RECO 1 
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PADDR  =  ADDR 
END  IF 
C 

C  BUILD  THE  CELL  MASK  FROM  THE  2  BIT  RECORDS  WITHIN  THE  LOGICAL 

C  RECORD;  00  -  BOTH  LAND  AND  WATER,  01  -  LAND,  10  -  WATER.  IF 

C  THE  2  BIT  RECORD  IS  00,  SET  THE  CELL  MASK  VALUE  TO  THE  LOGICAL 

C  ADDRESS  OF  THE  ONE  DEGREE  DATA  RECORD 

C 

SUM  =  0 
DO  20  I  =  0,9 
DO  10  J  =  0,9 

COMPUTE  THE  BIT  POSITION  WITHIN  THE  PHYSICAL  RECORD,  THE 
BYTE  POSITION  WITHIN  THE  PHYSICAL  RECORD,  AND  THE  BIT 
POSITION  WITHIN  THE  BYTE 

BPOS  =  (I*10+J)*2 
BYT  =  BPOS/a 
BIT  =  MOD (BPOS, 8) 

SET  =  MOD(ICHAR(IREC01(CELL+BYT) )/2»*BIT,4) 

IF  (SET.EQ.O)  THEN 
SUM  =  SUM+1 

CELMSK(J,I)  =  RECOl+SUM 
ELSE 

CELMSK(J,I)  =  SET 
END  IF 

10  CONTINUE 

20  CONTINUE 

PRECOl  =  RECOl 
END  IF 


GET  THE  ADDRESS,  OR  GET  THE  FLAG  IF  ALL  LAND  OR  WATER 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


LONP  «  MOD (LON, 10) 

LATP  «  MOD(LAT,10) 

IF  (CELMSK(LONP,LATP) .GT.2)  THEN 
ONEFLG  «  '  ' 

RECOO  «  CELMSK( LONP, LATP) 

ELSE 

ONEFLG  =  CTYPE (CELMSK( LONP, LATP ) ) 

END  IF 
RETURN 
END 

SUBROUTINE  ALLMSK  :  SET  THE  ENTIRE  MASK  TO  EITHER  ’L’  OR  'W'. 

AUTHOR  :  JAN  C.  DEPNER,  04/03/89 
ARGUMENTS  : 

MASK  -  (C)  ARRAY  CONTAINING  THE  ONE  MINUTE  GRID  OF  LAND 

MASK  POINTS  FOR  THE  ONE  DEGREE  CELL  THAT  CONTAINS 
THE  LATITUDE  AND  LONGITUDE  POINT  REQUESTED 
ALFLAG  -  (C)  ALL  LAND/WATER  FLAG  ('L'  OR  ’W) 

SUBROUTINE  ALLMSK ( ALFLAG , MASK ) 

IMPLICIT  INTEGER  (A-Z) 

CHARACTER*!  MASK(0 : 59 , 0 ; 59 ), ALFLAG 
DO  20  I  *  0,59 
DO  10  J  =  0,59 

MASK(J,I)  =  ALFLAG 
10  CONTINUE 
20  CONTINUE 
RETURN 
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END 


C 

C  SUBROUTINE  BLDMSK  ;  BUILDS  THE  MASK  FROM  THE  ONE  DEGREE  DATA 

C  RECORDS  IN  THE  DATA  BASE. 

C 

C  AUTHOR  :  JAN  C.  DEPNER,  04/03/89 

C 

C  ARGUMENTS  : 

C 

C  MASK 

C 
C 
C 

C  VARIABLES  : 

C 

C  PHYREC 

C 

C  IRECOO 

C 

C  CTYPE 

C  LPERP 

C  RECOO 

C 

C  RECOl 

C  PADDR 

C 

SUBROUTINE  BLDMSK (MASK) 

IMPLICIT  INTEGER  (A-Z) 

PARAMETER  (PHYR£C=1024) 

CHARACTER* 1  IRECOO ( 0 : PHYREC-l ) , MASK (0:59,0:59), CTYPE (0:1) 

COMMON  /ADDRS/  LPERP, RECOO, RECOl 
DATA  PADDR  /-!/  ,  CTYPE  /'W^’L'/ 

START  =  0 

COMPUTE  THE  FIRST  BYTE  POSITION  WITHIN  THE  PHYSICAL  RECORD,  AND  THE 
PHYSICAL  RECORD  ADDRESS 

POS  »  MOD (RECOO, LPERP) *32 
ADDR  *  INT(REC00/LPERP+1) 

IF  THE  ADDRESS  HAS  CHANGED,  READ  A  NEW  RECORD 

IF  ( ADDR. NE. PADDR)  THEN 
READ ( 10 , REC=ADDR ) IRECOO 
PADDR  -  ADDR 
END  IF 

MAIN  LOOP,  GETS  TYPE  AND  COUNT  FROM  BYTES.  FIRST  TWO  BITS  ARE  TYPE 
FLAGS;  00-1  BYTE  WATER,  01-1  BYTE  LAND,  10-2  BYTE  WATER, 
11-2  BYTE  LAND.  NEXT  6  OR  14  BITS  ARE  COUNT 


10  MBYTE  =  ICHAR ( IRECOO ( POS ) ) 

IF  (MBYTE. NE.O)  THEN 

IF  THIS  IS  A  TWO  BYTE  RECORD,  SET  DFLAG 

IF  ( MBYTE. GE, 128)  THEN 
MBYTE  =  MBYTE-128 
DFLAG  =  1 
ELSE 

DFLAG  =  0 
END  IF 
C 


(C)  ARRAY  CONTAINING  THE  ONE  MINUTE  GRID  OF  LAND 
MASK  POINTS  FOR  THE  ONE  DEGREE  CELL  THAT  CONTAINS 
THE  LATITUDE  AND  LONGITUDE  POINT  REQUESTED 


(I)  CONSTANT  SET  TO  NUMBER  OF  BYTES  IN  A  PHYSICAL 
RECORD  IN  THE  DATA  BASE. 

(C)  ARRAY  CONTAINING  AN  INPUT  RECORD  FROM  THE 
DATA  BASE. 

(C)  ARRAY  CONTAINING  THE  LAND  AND  WATER  FLAGS. 

(I)  LOGICAL  RECORDS  PER  PHYSICAL  RECORD. 

(I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  DATA  RECORD 
(NOT  USED  IN  THIS  ROUTINE). 

(I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  CELL  MAP. 
(I)  PREVIOUS  PHYSICAL  RECORD  ADDRESS. 
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GET  THE  TYPE  AND  COUNT 

TYPE  =  MBYTE/64 
COUNT  =  MBYTE-TYPE *64 

IF  THIS  IS  A  TWO  BYTE  RECORD,  GET  THE  REST  OF  THE  COUNT  FROM 
THE  NEXT  BYTE 

IF  (DFLAG.EQ.l)  THEN 

CALL  MOVPOS(POS, IRECOO,PADDR) 

COUNT  =  COUNT*256+ICHAR(IRECOO<POS) ) 

END  IF 

LOOP  THROUGH  THE  MASK  AND  SET  THE  CHARACTERS 

FINISH  =  START+COUNT-1 
DO  20  J  =  START, FINISH 
LATP  =  J/60 
LONP  =  MOD (J, 60) 

MASK ( LONP , LATP )  =  CTYPE(TYPE) 

20  CONTINUE 

MOVE  THE  BYTE  POSITION  AND  GO  TO  TOP  OF  LOOP 

CALL  MOVPOS(POS,IRECOO,PADDR) 

START  =  FINISH+1 
GO  TO  10 
END  IF 


FLIP  THE  TYPE  FLAG  AND  SET  THE  REMAINDER  OF  THE  MASK 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


25  TYPE  =  l-TYPE 

DO  30  J  =  START, 3599 
LATP  «  J/60 
LONP  «  MOD (J, 60) 

MASK (LONP, LATP)  =  CTYPE(TYPE) 
30  CONTINUE 

PTYPE  «  TYPE 

RETURN 

END 


SUBROUTINE  MOVPOS  ;  INCREMENTS  THE  BYTE  POSITION  AND  CHANGES 
TO  THE  NEXT  OVERFLOW  RECORD  IN  THE  CHAIN  IF  NECCESSARY. 

AUTHOR  .  JAN  C.  DEPNER,  04/03/89 


ARGUMENTS  : 


POS 

IRECOO 

PADDR 


(I)  BYTE  POSITION  WITHIN  THE  PHYSICAL  RECORD. 
(C)  ARRAY  CONTAINING  AN  INPUT  RECORD  FROM  THE 
DATA  BASE. 

(I)  PREVIOUS  PHYSICAL  RECORD  ADDRESS. 


VARIABLES  ; 


PHYREC 

LPERP 

RECOO 

RECOl 


(I)  CONSTANT  SET  TO  NUMBER  OF  BYTES  IN  A  PHYSICAL 
RECORD  IN  THE  DATA  BASE. 

(1)  LOGICAL  RECORDS  PER  PHYSICAL  RECORD. 

(I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  DATA  RECORD 
(NOT  USED  IN  THIS  ROUTINE). 

(I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  CELL  MAP. 


SUBROUTINE  MOVPOS(POS, IRECOO, PADDR) 
IMPLICIT  INTEGER  (A-Z) 
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PARAMETER  ( PHyREC=1024 ) 

CHARACTER*!  IRECOO ( 0 ; PHYREC-1 ) 

COMMON  /ADDRS/  LFERP , RECOO , RECOl 
POS  =  POS+1 

IF  THIS  IS  THE  END  OF  THE  LOGICAL  RECORD,  GET  THE  OVERFLOW  POINTER 
FROM  THE  LAST  2  BYTES  OF  THE  RECORD 

IF  (MOD(POS+2,32) .EQ.O)  THEN 

LREC  =  ICHAR(IRECOO(POS) )*2S6+ICHAR(IRECOO(POS+1) )+REC01 

IF  LREC  =  RECOl,  THERE  IS  NO  OVERFLOW  POINTER  AND  THIS  IS  THE 
END  OF  THE  CHAIN 

IF  ( LREC. NE. RECOl)  THEN 

COMPUTE  THE  PHYSICAL  ADDRESS,  IF  IT  HAS  CHANGED,  READ  A  NEW 
RECORD 

ADDR  =  LREC/LPERP+1 
IF  (ADDR.NE.PADDR)  THEN 
READ(10,REC=ADDR) IRECOO 
PADDR  =ADDR 
END  IF 

COMPUTE  THE  NEW  BYTE  POSITION  WITHIN  THE  PHYSICAL  RECORD 

POS  =  MOD(LREC,LPERP)*32 
END  IF 
END  IF 
RETURN 
END 
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APPENDIX  B 


B-1 


/include  <math.h> 

/include  <fcntl.h> 

/include  <syB\types. h> 

/include  <sys\stat.h> 

/include  <stdio.h> 

/include  <io.h> 

/define  PHYREC  1024 
/define  LPERP  PHYREC/32 
/define  FILENAME  "landmaak.dat" 

/define  SIGN_OF(x)  ((x)<0.0  ?  -1  :  1) 

/* 

FUNCTION  MASKOl  :  RETURNS  A  LOGICAL  FLAG  INDICATING  THAT  THE 
ONE  MINUTE  GRID  NEAREST  TO  THE  LATITUDE  AND  LONGITUDE  POSITION 
PASSED  TO  THE  SUBROUTINE  IS  EITHER  ON  LAND  OR  WATER.  THE  POINT 
MAY  BE  DEFINED  AS  DEGREES  DECIMAL  OR  DEGREES  AND  MINUTES  DECIMAL. 
THE  FLAG  RETURNED  IS  SET  TO  .TRUE.  FOR  LAND  OR  .FALSE.  FOR  WATER. 
THIS  ROUTINE  CALLS  MASK60  TO  GET  THE  DEGREE  SQUARE  MrtSK. 

THIS  PROGRAM  USES  A  DATA  BASE  DESIGNED  BY: 


JAN  C.  DEPNER 
NAVAL  OCEANOGRAPHIC  OFFICE 
SYSTEMS  TECHNOLOGY  BRANCH 
STENNIS  SPACE  CENTER,  MS 
39522-5001 


JAMES  A.  HAMMACK 

NAVAL  OCEAN  R&D  ACTIVITY 

CODE  117 

STENNIS  SPACE  CENTER,  MS 
39529- 


AUTHOR  :  JAN  C.  DEPNER,  04/03/89 


ARGUMENTS  : 

latd  -  (F)  LATITUDE  DEGREES  OF  THE  POINT  (SOUTH  NEGATIVE) 

latm  -  (F)  LATITUDE  MINUTES  OF  THE  POINT  (UNSIGNED  UNLESS 

LATD  IS  0.0) 

lond  -  (F)  LONGITUDE  DEGREES  OF  THE  POINT  (WEST  NEGATIVE  OR 

0  TO  360) 

lonra  -  (F)  LONGITUDE  MINUTES  OF  THE  POINT  (UNSIGNED  UNLESS 

LOND  IS  0.0) 

VARIABLES  : 


mask 


latp 

lonp 

lat 

Ion 


-  (C)  ARRAY  CONTAINING  THE  ONE  MINUTE  GRID  OF  LAND 
MASK  POINTS  FOR  THE  ONE  DEGREE  CELL  THAT  CONTAINS 
THE  LATITUDE  AND  LONGITUDE  POINT  REQUESTED 

-  (I)  LATITUDE  INDEX  INTO  MASK  ARRAY 

-  (I)  LONGITUDE  INDEX  INTO  MASK  ARRAY 

-  (F)  LATITUDE  IN  DECIMAL  DEGREES 

-  (F)  LONGITUDE  IN  DECIMAL  DEGREES 


char  maskOl ( latd, latm, lond, lonm) 
float  latd,  latm,  lond,  lonm; 

{ 

static  unsigned  char  mask [3600]; 
static  char  alflag,  flag60  =  127; 
float  lat.  Ion; 
double  dummy; 

int  dlat,  dlon,  latp,  lonp; 
static  int  plat  =  -1,  plon  =  -1; 

void  ma8k60 ( ) ; 
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/* 


*/ 


/* 


*/ 


/* 


*/ 


/* 


*/ 


/* 


*/ 


/* 


V 


> 

/* 


CONVERT  LAT  AND  LON  TO  DEGREES  DECIMAL,  ADJUST  TO  NEAREST  MINUTE 

lat  =  latd  +  (latm/60.0)  •  SIGN_OF( latd)  +  90,0; 

Ion  =  lond  +  (lonm/60.0)  •  SIGN_OF( lend) ; 
if  (lon<0.0)  Ion  +=  360.0; 
lat  +=  .00833333333333; 

Ion  +=  .00833333333333; 

IF  THIS  IS  IN  THE  SAME  ONE  DEGREE  SQUARE  AS  THE  LAST  ACCESS, 

DO  NOT  READ  A  NEW  MASK 

dlat  =  (int)lat; 
dlon  =  (int)lon; 
if  (dlon!=pion  ||  dlat!=plat)  { 

GET  THE  ONE  DEGREE  MASK 

maBk60(mask,  lat-90.0.  Ion, Sif  lag60) ; 

} 

IF  THE  ALL  LAND  OR  WATER  FLAG  IS  NOT  SET,  GET  THE  SINGLE  POINT 
FROM  THE  MASK 

if  {flag60<0)  { 

latp  =  modf ( (double) lat, fidummy) *60,0; 
lonp  =  modf ( (double) Ion, fidummy) *60  0; 
alflag  =  maBk( ( int ) latp*60+ ( int ) lonp) ; 

} 

elBe  { 

alflag  =  flag60; 

> 

SAVE  THE  PREVIOUS  LAT  AND  LON  DEGREES 

plat  =  dlat; 
plon  =  dlon; 

SET  THE  LOGICAL  VALUE  BASED  ON  THE  SETTING  OF  ALFLAG 

if  (alflag)  { 
return  ( 1 ) ; 

} 

elBe  { 

return  (0) ; 

} 


FUNCTION  MASK60  :  GIVEN  A  REAL  LATITUDE  AND  LONGITUDE  IN  DEGREES 
(WEST  LONGITUDE  NEGATIVE  OR  0  TO  360  LONGITUDE  AND  SOUTH  LATITUDE 
NEGATIVE)  THIS  SUBROUTINE  RETURNS  A  CHARACTER  ARRAY  CONTAINING  A 
ONE  MINUTE  LAND  MASK  FOR  THE  ONE  DEGREE  CELL  THAT  THE  LATITUDE 
AND  LONGITUDE  POINT  IS  IN.  THE  CHARACTER  ARRAY  CONTAINS  3600 
CHARACTERS  AND  IS  SET  UP  AS  'UNSIGNED  CHT^ ( 60 ] [ 60 ] ’ .  EACH 
CHARACTER  IS  SET  TO  EITHER  ' 1 '  OR  ' 0 ' .  IF  ALL  THE  POINTS  ARE 
LAND  OR  ALL  ARE  WATER,  THE  1  CHARACTER  FLAG  'ALFLAG'  WILL  BE  SET 
TO  THE  CORRESPONDING  VALUE.  IF  BOTH  LAND  AND  WATER  POINTS  ARE 
PRESENT,  'ALFLAG'  WILL  BE  SET  TO  '-1'.  THE  MASK  IS  SET  UP  WITH 
POINT  'MASK(0)(0)'  BEING  THE  SOUTHWEST  CORNER  OF  THE  CELL  AND  POINT 
■MASK(59] [59] '  THE  NORTHEAST  CORNER.  THE  FIRST  SUBSCRIPT  REPRESENTS 
THE  LONGITUDE  MINUTES,  THE  SECOND  IS  LATITUDE  MINUTES  (WITH  THE 
COORDINATE  SYSTEM  SET  UP  AS  0  -  360  LONGITUDE  AND  0  -  180  LATITUDE). 
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THIS  PROGRAM  USES  A  DATA  BASE  DESIGNED  BY: 


JAN  C.  DEPNER 
NAVAL  OCEANOGRAPHIC  OFFICE 
SYSTEMS  TECHNOLOGY  BRANCH 
STENNIS  S6ACE  CENTER,  MS 
39522-5001 


JAMES  A.  HAMMACK 

NAVAL  OCEAN  R&D  ACTIVITY 

CODE  117 

STENNIS  SPACE  CENTER,  MS 
39529- 


AUTHOR  ;  JAN  C.  DEPNER,  04/03/89 
ARGUMENTS  : 

mask  -  (C)  ARRAY  CONTAINING  THE  ONE  MINUTE  GRID  OF  LAND 

MASK  POINTS  FOR  THE  ONE  DEGREE  CELL  THAT  CONTAINS 
THE  LATITUDE  AND  LONGITUDE  POINT  REQUESTED 
lat  -  (F)  LATITUDE  DEGREES  OF  THE  POINT  (0  -  180) 

Ion  -  (F)  LONGITUDE  DEGREES  OF  THE  POINT  (0  -  360) 

alflag  -  (C)  ALL  LAND/WATER  FLAG  RETURNED 

VARIABLES  : 


*/ 


PHYREC 

ireclO 

LPEkP 

recOO 

recOl 

paddr 


(I)  CONSTANT  SET  TO  NUMBER  OF  BYTES  IN  A  PHYSICAL 
RECORD  IN  THE  DATA  BASE. 

-  (C)  ARRAY  CONTAINING  AN  INPUT  RECORD  FROM  THE 
DATA  BASE. 

-  (I)  LOGICAL  RECORDS  PER  PHYSICAL  RECORD. 

-  (I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  DATA  RECORD 
(NOT  USED  IN  THIS  ROUTINE). 

-  (I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  CELL  MAP. 

-  (I)  PREVIOUS  PHYSICAL  RECORD  ADDRESS. 


void  ma8k60(ma8k, lat, Ion, alflag) 
char  ma8k[],  *alflag; 
float  lat.  Ion; 

{ 

static  long  recOO,  recOl; 

static  int  first  =  1,  plat  =  -1,  plon  =  -1,  lunfil; 
int  clat,  cion; 

void  rdten(),  rdone(),  allmskO,  bldmsk(); 

/* 

IF  FIRST  TIME  THROUGH,  OPEN  FILE  AND  SET  LPERP 

*/ 

if  (first)  { 

lunfil  =  open ( FILENAME, 0_BINARy jo_RDONLY) ; 
first  =  0; 

} 

/* 

IF  LAT  >=  90  SET  TO  ALL  WATER  AND  RETURN 

*/ 

if  (lat  >=  90.0)  { 

*alflag  =  0; 

allm8k( *alf lag, mask) ; 

return; 

} 

/* 

CALCULATE  SOUTHWEST  CORNER  OF  CELL 

*/ 

clat  «  (int) (lat+90.0) ; 
if  (lon<0.0)  { 

cion  =  (int) (lon+360.0) ; 

} 


B-4 


/* 


*/ 


/* 


*/ 


/* 


*/ 


/* 


*/ 


/* 


*/ 


} 

/* 


else  { 

cion  =  ( ( int ) Ion) %360; 

} 


IF  THE  LAT  AND  LON  ARE  NOT  IN  THE  SAME  ONE  DEGREE  SQUARE  AS  THE 
PREVIOUS  LAT  AND  LON,  OR  ALFLAG  =  127  (FIRST  CALL  FROM  MASKOl), 
PROCESS  THE  CELL 

if  (clonl=plon  ||  clatl=plat  ||  *alf lag==127 )  { 

*alflag  =  '  ' ; 

READ  THE  TEN  DEGREE  SQUARE  MAP  AND  CHECK  FOR  ALL  LAND  OR  WATER 

rdten( lunfil , &rec01 , clat , cion, al flag) ; 
if  (*alflag>=0)  { 

allmsk( *alf lag, mask) ; 

} 

else  { 

READ  THE  ONE  DEGREE  CELL  MAP  AND  CHECK  FOR  ALL  LAND  OR  WATER 

rdone ( lunfil , recOl , firecOO, clat , cion, alf lag ) ; 
if  (*alflag>=0)  { 
allmsk( *alf lag, mask) ; 

} 

else  { 

READ  THE  RUN  LENGTH  DATA  FOR  THE  ONE  DEGREE  CELL  AND  BUILD 
THE  MASK 

bldroek(lunfil,rec00,rec01,mask) ; 

} 


SAVE  THE  LAT  AND  LON  DEGREES 

plat  =  clat; 
plon  cion; 


FUNCTION  RDTEN  :  READS  THE  TEN  DEGREE  SQUARE  MAP  AT  THE 
BEGINNING  OF  THE  DATA  BASE.  RETURNS  A  FLAG  INDICATING  IF  THE 
ENTIRE  SQUARE  WAS  EITHER  LAND  OR  WATER.  IF  NOT,  THE  VARIABLE 
•RECOl’  IS  SET  TO  THE  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  CELL  MAP. 
THE  FLAG  RETURNED  IS  SET  TO  ’1’  FOR  LAND  OR  '0'  FOR  WATER. 


AUTHOR  ;  JAN  C.  DEPNER,  04/03/89 
ARGUMENTS  : 


lunfil 

-  (I) 

recOl 

-  (L) 

lat 

-  (F) 

Ion 

-  (F) 

tenf Ig 

-  (C) 

VARIABLES 

; 

PHYREC 

-  (I)  ' 

RECORD 

ireclO 

-  (C) 

DATA  B. 

INPUT  FILE  HANDLE 

RECORD  NUMBER  OF  ONE  DEGREE  CELL  MAP 
LATITUDE  DEGREES  OF  THE  POINT  (0  -  180) 
LONGITUDE  DEGREES  OF  THE  POINT  (0  -  360) 
ALL  LAND /WATER  FLAG  RETURNED 
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LPERP  -  (I)  LCXJICAL  RECORDS  PER  PHYSICAL  RECORD. 

recOO  -  (I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  DATA  RECORD 

(NOT  USED  IN  THIS  ROUTINE). 

recOl  -  (I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  CELL  MAP. 

paddr  -  (I)  PREVIOUS  PHYSICAL  RECORD  ADDRESS. 

*/ 

void  rdten(lunf il,rec01, lat, lon,tenflg) 
int  lunfil,  lat,  Ion; 
long  *rec01; 
char  *tenflg; 

{ 

static  long  paddr  =  -1,  pword  =  -1; 
static  unsigned  char  ireclO[PHYREC] ; 
long  addr,  word,  Irec,  pnt,  Istat; 

/* 

COMPUTE  THE  WORD  (4  BYTE)  POSITION  WITHIN  THE  648  WORD  TEN 
DEGREE  SQUARE  MAP,  MAP  STARTS  IN  THE  2ND  LOGICAL  RECORD 

*! 

word  =  (lat/10)*36+lon/10; 

/* 

IF  THE  WORD  POSITION  HAS  CHANGED  GET  THE  NEW  ONE  DEGREE  CELL 
LOGICAL  ADDRESS 

*/ 

if  (word i=pword)  { 

/* 

COMPUTE  THE  LOGICAL  AND  PHYSIC7VL  ADDRESS  OF  THE  ONE  DEGREE 
CELL  POINTER  WITHIN  THE  TEN  DEGREE  MAP 

*/ 

Irec  =  (word*4) /32+1; 

addr  =  ( Irec/ (long) (LPERP) )*PHYREC; 

/* 

IF  THE  PHYSICAL  ADDRESS  HAS  CHANGED,  READ  A  NEW  RECORD 

*/ 

if  (addr i=paddr)  { 

Istat  =  1 seek (lunfil, addr, SEEK_SET ) ; 

Istat  =  read ( lunfil, ireclO,PHYR£C) ; 
paddr  =  addr; 

} 

f* 

COMPUTE  THE  BYTE  POSITION  OF  THE  ONE  DEGREE  CELL  LOGICAL 
ADDRESS  WITHIN  THE  PHYSICAL  RECORD 

*/ 

pnt  =  ( (word+8)»4)%PHYREC; 

/* 

BUILD  THE  LOGICAL  ADDRESS  FROM  THE  4  BYTES  AT  'PNT' 

*! 

*rec01  =  irecl0lpnt]*16777216+irecl0lpnt+l)*65536+ 
irecl0(pnt+2 ]*256+irecl0(pnt+3] ; 
pword  =  word; 

} 

/* 

IF  THE  TEN  DEGREE  SQUARE  IS  ALL  LAND  OR  WATER,  SET  THE  FLAG 

*/ 

if  (*rec01<3)  { 

•tenflg  =  (int)*rec01%2; 

} 

else  { 

♦tenflg  =  -1; 

} 

return  ; 
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} 

/* 

FUNCTION  RDONE  :  READS  THE  ONE  DEGREE  CELL  MAP.  RETURNS  A 
FLAG  THAT  INDICATES  IF  THE  ENTIRE  SQUARE  IS  EITHER  LAND  OR 
WATER.  IF  NOT,  THE  VARIABLE  'RECOO*  IS  SET  TO  THE  LOGICAL 
ADDRESS  OF  THE  ONE  DEGREE  DATA  RECORD.  THE  FLAG  RETURNED  IS 
SET  TO  1  FOR  LAND  OR  0  FOR  WATER. 

AUTHOR  :  JAN  C.  DEPNER,  04/03/89 

ARGUMENTS  : 

lat  -  (R)  LATITUDE  DEGREES  OF  THE  POINT  (0  -  180) 

Ion  -  (R)  LONGITUDE  DEGREES  OF  THE  POINT  (0  -  360) 

oneflg  -  (C)  ALL  LAND/WATER  FLAG  RETURNED 

VARIABLES  : 

PHYREC  -  (I)  CONSTANT  SET  TO  NUMBER  OF  BYTES  IN  A  PHYSICAL 

RECORD  IN  THE  DATA  BASE. 

celmsk  -  (I)  ARRAY  CONTAINING  THE  LOGICAL  ADDRESSES  OF 

THE  ONE  DEGREE  CELL  DATA  RECORDS. 
irecOl  -  (C)  ARRAY  CONTAINING  AN  INPUT  RECORD  FROM  THE 

DATA  BASE. 

LPERP  -  (I)  LOGICAL  RECORDS  PER  PHYSICAL  RECORD. 

recOO  -  (I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  DATA  RECORD 

(NOT  USED  IN  THIS  ROUTINE). 

recOl  -  (I)  LOGICAL  ADDRESS  OF  T»iE  ONE  DEGREE  CELL  MAP. 

paddr  -  (I)  PREVIOUS  PHYSICAL  Ri:CORD  ADDRESS. 

*/ 

void  rdone (lunf 11, rec01,rec00, lat, Ion, oneflg) 
int  lunfil,  lat.  Ion; 
long  recOl,  *rec00; 
char  *oneflg; 

{ 

static  long  celmsk ( 10] ( 10] ,  paddr  =  -1,  precOl  =  -1; 
static  unsigned  char  irecOl [PHYREC] ; 
static  int  power []  =  {1,2,4,8,16,32,64,128}; 
long  addr,  cell,  Istat; 

int  sum,  bpos,  byt,  bit,  set,  i,  j,  latp,  lonp; 

/* 

COMPUTE  THE  CELL  POSITION  WITHIN  THE  PHYSICAL  RECORD  AND  THE 
PHYSICAL  RECORD  ADDRESS 

*/ 

cell  =  (rec01%(long) (LPERP) )*32; 
addr  =  (rec01/(long) (LPERP) )*PHYREC; 

/* 

IF  LOGICAL  RECORD  HAS  CHANGED,  BUILD  A  NEW  CELL  MASK 

*/ 

if  ( recOl l=prec01)  { 

/* 

IF  THE  ADDRESS  HAS  CHANGED,  READ  A  NEW  RECORD 

*/ 

if  (addr 1 -paddr)  { 

istat  =  lseek(iunfil,addr,0) ; 

Istat  »  read( iunfil, irecoi, PHYREC) ; 
paddr  *  addr; 

} 

I* 

BUILD  THE  CELL  MASK  FROM  THE  2  BIT  RECORDS  WITHIN  THE  LOGICAL 
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RECORD;  00  -  BOTH  LAND  AND  WATER,  01  -  LAND,  10  -  WATER.  IF 
THE  2  BIT  RECORD  IS  00,  SET  THE  CELL  MASK  VALUE  TO  THE  LOGICAL 
ADDRESS  OF  THE  ONE  DEGREE  DATA  RECORD 

*/ 

sum  =0; 

for  (i=0;  i<=9;  i++)  { 
for  {j=0;  j<=9;  j++)  { 

/* 

COMPUTE  THE  BIT  POSITION  WITHIN  THE  PHYSICAL  RECORD,  THE 
BYTE  POSITION  WITHIN  THE  PHYSICAL  RECORD,  AND  THE  BIT 
POSITION  WITHIN  THE  BYTE 

*/ 

bpos  =  (i*10+j)*2; 
byt  =  bpos/8; 
bit  =  bpoB%8; 

set  =  ( irec01(cell+byt] /power [ bit ) )%4; 
if  (Iset)  { 

Bum++; 

celmsk[j)(i]  =  recOl+sum; 

} 

else  { 

celmBk(j](i]  =  set; 

} 

} 

} 

precOl  =  recOl; 

} 

/* 

GET  THE  ADDRESS,  OR  GET  THE  FLAG  IF  ALL  LAND  OR  WATER 

*/ 

lonp  =  Ion%10; 
latp  »  lat%10; 
if  (celm8k[ lonp] [ latp]>2 )  < 

•oneflg  =  -1; 

*rec00  »  celmsk] lonp] ( latp] ; 

} 

else  { 

*oneflg  -  celmsk] lonp] [ latp}%2; 

} 

return; 

} 

/* 

FUNCTION  ALLMSK  ;  SET  THE  ENTIRE  MASK  TO  EITHER  ‘L’  OR  ’W. 

AUTHOR  :  JAN  C.  DEPNER,  04/03/89 
ARGUMENTS  : 

mask  -  (C)  ARRAY  CONTAINING  THE  ONE  MINUTE  GRID  OF  LAND 

MASK  POINTS  FOR  THE  ONE  DEGREE  CELL  THAT  CONTAINS 
THE  LATITUDE  AND  LONGITUDE  POINT  REQUESTED 
alflag  -  (C)  ALL  LAND/WATER  FLAG 

*/ 

void  allmsk(alf lag, mask) 
char  alflag,  mask]]; 

{ 


int  i ,  j ; 

for  (i=0;  i<=59;  i++)  { 
for  (j=0;  j<*59;  j++)  { 
mask[i*60+j]  =  alflag; 


B-8 


} 

> 

return; 

} 

/* 

FUNCTION  BLDMSK  ;  BUILDS  THE  MASK  FROM  THE  ONE  DEGREE  DATA 
RECORDS  IN  THE  DATA  BASE. 

AUTHOR  :  JAN  C.  DEPNER,  04/03/89 

ARGUMENTS  ; 

mask  -  (C)  ARRAY  CONTAINING  THE  ONE  MINUTE  GRID  OF  LAND 

MASK  POINTS  FOR  THE  ONE  DEGREE  CELL  THAT  CONTAINS 
THE  LATITUDE  AND  LONGITUDE  POINT  REQUESTED 

VARIABLES  : 

PHYREC  -  (I)  CONSTANT  SET  TO  NUMBER  OF  BYTES  IN  A  PHYSICAL 

RECORD  IN  THE  DATA  BASE. 

irecOO  -  (C)  ARRAY  CONTAINING  AN  INPUT  RECORD  FROM  THE 

DATA  BASE. 

LPERP  -  (I)  LOGICAL  RECORDS  PER  PHYSICAL  RECORD. 

recOO  -  (I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  DATA  RECORD 

(NOT  USED  IN  THIS  ROUTINE). 

recOl  -  (I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  CELL  MAP. 

paddr  -  (I)  PREVIOUS  PHYSICAL  RECORD  ADDRESS. 

*/ 

void  bldmski lunf il,rec00, recOl^mask) 

Int  lunfil; 
long  recOO,  recOl; 
char  ma8k[ ] ; 

{ 

static  long  paddr  »  -1; 
static  unsigned  char  irecOO [ PHYREC ] ; 
int  start,  type,  count,  finish,  j; 
long  addr,  pos,  Istat; 
unsigned  char  mbyte,  dflag; 

void  movpos ( ) ; 

start  «  0; 

/* 

COMPUTE  THE  FIRST  BYTE  POSITION  WITHIN  THE  PHYSICAL  RECORD,  AND  THE 
PHYSICAL  RECORD  ADDRESS 

*/ 

pos  =  (rec00% ( long) (LPERP) ) *32 ; 
addr  =  (recOO/ ( long) (LPERP) ) *PHYREC; 

/* 

IF  THE  ADDRESS  HAS  CHANGED,  READ  A  NEW  RECORD 

*! 

i f  ( addr 1 »paddr )  { 

printf  ("paddr  *  %d  %o,  addr  =  %d  %o\n",paddr,paddr, addr, addr) ; 

Istat  «  l8eek(lunfil,addr,0) ; 

Istat  »  read ( lunfil , irecOO, PHYREC) ; 
paddr  >  addr; 

} 

I* 

MAIN  LOOP,  GETS  TYPE  AND  COUNT  FROM  BYTES.  FIRST  TWO  BITS  ARE  TYPE 
FLAGS;  00-1  BYTE  WATER,  01-1  BYTE  LAND,  10-2  BYTE  WATER, 
11-2  BYTE  LAND.  NEXT  6  OR  14  BITS  ARE  COUNT 
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mbyte  =  irecOO[pos); 
while  (mbyte)  { 

IF  THIS  IS  A  TWO  BYTE  RECORD,  SET  DFLAG 

if  (mbyte>=128)  { 

mbyte  =  mbyte- 12 8; 
dflag  =  1; 

} 

else  { 

dflag  =  0; 

} 

GET  THE  TYPE  AND  COUNT 

type  =  mbyte/64; 
count  =  mbyte-type*64; 

IF  THIS  IS  A  TWO  BYTE  RECORD,  GET  THE  REST  OF  THE  COUNT  FROM 
THE  NEXT  BYTE 

if  (dflag)  { 

movpo8(  lunf  il,  SipoB,  irec00,&paddr,  recOl )  ; 
count  =  count*256+irec00(poa ) ; 

} 

LOOP  THROUGH  THE  MASK  AND  SET  THE  CHARACTERS 

finish  =  start+count-l ; 
for  (j=start;  j<=fini8h;  j++)  { 
ma8k[j]  =  type; 

} 

MOVE  THE  BYTE  POSITION  AND  GO  TO  TOP  OF  LOOP 

movpo8(lunfil,&po8, irecOO, &paddr, recOl ) ; 
start  »  finish+l; 
mbyte  -  irec00[po8]; 


FLIP  THE  TYPE  FLAG  AND  SET  THE  REMAINDER  OF  THE  MASK 
type  =  1-type; 

for  (j=8tart;  j<=3599;  j++)  { 
maek[j]  =  type; 

} 

return; 


FUNCTION  MOVPOS  :  INCREMENTS  THE  BYTE  POSITION  AND  CHANGES 
TO  THE  NEXT  OVERFLOW  RECORD  IN  THE  CHAIN  IF  NECCESSARY. 

AUTHOR  :  JAN  C.  DEPNER,  04/03/89 

ARGUMENTS  : 

pos  -  (I)  BYTE  POSITION  WITHIN  THE  PHYSICAL  RECORD. 

irecOO  -  (C)  ARRAY  CONTAINING  AN  INPUT  RECORD  FROM  THE 

DATA  BASE. 

paddr  -  (I)  PREVIOUS  PHYSICAL  RECORD  ADDRESS. 


VARIABLES 


PHYREC 


-  (I)  CONSTANT  SET  TO  NUMBER  OF  BYTES  IN  A  PHYSICAL 

RECORD  IN  THE  DATA  BASE. 

LPERP  -  (I)  LOGICAL  RECORDS  PER  PHYSICAL  RECORD. 

recOO  -  (I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  DATA  RECORD 

(NOT  USED  IN  THIS  ROUTINE). 

recOl  -  (I)  LOGICAL  ADDRESS  OF  THE  ONE  DEGREE  CELL  MAP. 

*/ 

void  movpo8(lunfil,po8, irec00,paddr,rec01) 
int  lunfil; 

un8igned  char  irecOO[]; 
long  *poa,  *paddr,  recOl; 

{ 

long  Irec,  addr,  latat; 

(*po8)++; 

/* 

IF  THIS  IS  THE  END  OF  THE  LOGICAL  RECORD,  GET  THE  OVERFLOW  POINTER 
FROM  THE  LAST  2  BYTES  OF  THE  RECORD 

V 

if  ( (*po8+2)%32==0)  { 

Irec  =  irecOO( *po8] *256+irec00( *poB+l )+rec01; 

/* 

IF  LREC  =  RECOl,  THERE  IS  NO  OVERFLOW  POINTER  AND  THIS  IS  THE 
END  OF  THE  CHAIN 

*/ 

if  ( Irec 1 =rec01 )  { 

t* 

COMPUTE  THE  PHYSICAL  ADDRESS,  IF  IT  HAS  CHANGED,  READ  A  NEW 
RECORD 

*/ 

addr  =  ( Irec/ ( long) (LPERP) ) ‘PHYREC; 
if  (addr l=*paddr )  { 
latat  >  l8eelc(lunfil,addr,0) ; 
l8tat  -  read(lunfil,irecOO,PHYREC) ; 

*paddr  >  addr; 

} 

/• 

COMPUTE  THE  NEW  BYTE  POSITION  WITHIN  THE  PHYSICAL  RECORD 

*/ 

*po8  =  ( lrec%( long) (LPERP) ) *32; 

} 

} 

return; 

} 
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